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Abstract Equation of motion for a comet in the Oort cloud is numerically solved. 
Orbital evolution of the comet under the action of the gravity of the Sun and the 
Galaxy is presented for various initial conditions. 

Oscillations of the Sun with respect to the galactic equatorial plane are taken 
into account. Real values of physical quantities concerning the gravitational action of 
the galactic neighbourhood of the Sun are important. The results are compared with 
currently used more simple models of the galactic tide. It turns out that physically 
improved models yield results which significantly differ from the results obtained on 
the basis of the conventional models. E.g., the number of returns of the comets into the 
inner part of the Solar System are about two times greater than it is in the conventional 
models. 

It seems that a comet from the Oort cloud can be a source of the dinosaurs extinc- 
tion at about 65 Myr ago. A close encounter of a star or an interstellar cloud disturbed 
a comet of the Oort cloud in the way that its semi-major axis increased/decreased 
above the value 5 x 10 4 AU and the comet hit the Earth. 

Keywords Oort cloud • Galaxy ■ Comets ■ Equation of motion ■ Orbital evolution 
1 Introduction 

Global galactic gravitational field influences motion of a comet in the Opik-Oort cloud 
(Opik 1932, Oort 1950) in the form of the galactic tide. The motion of the comet with 
respect to the Sun is important in better understanding of the Oort cloud. This paper 
deals with detailed numerical solution of the equation of motion of the comet under 
the gravity of the Sun and the galactic tide. We consider equation of motion presented 
by Klacka (2009a, 2009b). The x— and y— components of the acceleration come not 
only from the x— and y— components of the position of the comet, but also from the 
z— component of the position due to the gravity of the galactic disk. Our paper provides 
detailed numerical solutions of the equation of motion. The importance of the /"-terms 
is tested. Comparison with the conventional models of galactic tide is discussed (see, 
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e.g., Mihalas 1968, Heilser and Tremaine 1986, Levison et al. 2001, Dybczynski et al. 
2008). 



2 Orbital elements 

We discuss several models of galactic tide. Each of them is represented by an equation 
of motion. The equation of motion is numerically solved. On the basis of the positional 
and velocity vectors at each position, we can find orbital elements and their evolution. 
We can use Eq. (47) in Klacka (2004) for this purpose (the right-hand side of the last 
equation in Eq. 47 contains 1/e instead of 1). We can summarize the equations in 
the following form [osculating orbital elements: a - semi-major axis; e - eccentricity; 
i - inclination of the orbital plane to the reference plane - galactic equatorial plane; 
Q - longitude of the ascending node; cu - longitude of pericenter (the argument of 
pericenter/perihelion); q - perihelion distance; Q - aphelion distance; is the position 
angle of the particle on the orbit, when measured from the ascending node in the 
direction of the particle's motion, = u> + /]: 
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3 Current models 

3.1 Simple model 

The most simple model considering galactic tide is described by the well-known equa- 
tion of motion (see, e. g., Mihalas 1968, Heisler and Tremaine 1986) 
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C - [AtvG g + 2 (A 2 - B 2 )] C 
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where G is the gravitational constant, M© is mass of the Sun, p is the mass density of 
the Galaxy and A, B are Oort constants. Numerical solution of the equation of motion 
was presented by Pretka and Dybczynski (1994) and secular evolution by Klacka and 
Gajdosfk (2001): 
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3.2 Standard model 



Previous section presented the most simple action of galactic tide. A more elaborated 
model is often used (Heisler and Tremaine 1986, Levison et al. 2001, Dybczynski et al. 



4 



2008): 



ii 
dt 2 



d 2 v GM e 



— 3^ £ + {A-B)[A + B + 2A cos (2 u t)] £ 
2A(A - B) sin (2 u t) r) 



■q - 2A(A - B) sin (2 w *) £ 



dt 2 r3 

+ (A - B) [A + B - 2A cos (2 u t)] rj 

5 = -^C~ \_A-k G g + 2 (A 2 — -B 2 )] C 

r = \/£ 2 + ^ + C 2 , 
u; = A - B , (4) 

where G is the gravitational constant, M© is the mass of the Sun and the numerical 
values of the other relevant quantities are 

A — 13.0 km s 1 kpc 1 , 
B — — 13.0 km s _1 kpc 1 , 

g = 0.10 M© pc~ 3 . (5) 




t [10 years] t [10 years] 

Fig. 1 Orbital evolution of the comet situated in the Oort cloud with Oj n = 5 X 10 4 AU, 
ei n ra 0, iin = 90° under the influence of the solar gravity and the galactic tide based on the 
standard model for various mass densities in the neighbourhood of the Sun: p = 0.075 Mq 
pc~ 3 (dotted line), p = 0.1 M pc -3 (dashed line), p = 0.15 Mq pc~ 3 (solid line). 
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Fig. 1 represents the numerical solution of Eqs. (4) and (5). It shows the orbital 
evolution of the comet situated in the Oort cloud of comets (dj„ = 5 x 10 4 AU) with 
initial almost circular orbit perpendicular to the plane of the galactic equator (e in ?a 
0, ii n = 90°). Solar gravity and the galactic tide in terms of the standard model are 
considered for various local densities p in the neighbourhood of the Sun. 

In the literature the values of the local density span from p — 0.075 Mq pc~ 3 
(Creze et al. 1998) to p = 0.185 M pc~ 3 (Bahcall 1984). As it is shown in Fig. 1, the 
orbital evolution of the comet initially situated in the Oort cloud is very sensitive to 
the local density. For comparison we used the values p = 0.075 Mq pc -3 , p = 0.1 Mq 
pc~ 3 (Lcvison et al. 2001) and p = 0.15 M pc~ 3 (Wiegert and Tremaine 1999). 



3.3 Comparison of the current models 




t[10 8 years] 1 1 1 s years | 

Fig. 2 Orbital evolution of the comet situated in the Oort cloud with cii n = 5 X 10 4 AU, 
ei n « 0, ii n = 90° under the influence of the solar gravity and the galactic tide for the simple 
(dotted line) and the standard (solid line) models. Local density p = 0.1 Mq pc~ 3 is considered 
for both models. 



Figs. 2 and 3 compare orbital evolution of the comet calculated from the models 
discussed above. 

Fig. 2 shows the difference between the simple and the standard model for the the 
comet with initial values a,i n = 5 x 10 4 AU, ei n w and ii n = 90°. For the standard 
model (solid line) the number of returns of the comet to the inner part of the Solar 
System (perihelion distance less than 150 AU, approximately) is 5, while for the simple 
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t [10 s years] t [10 s years] 

Fig. 3 Orbital evolution of the comet situated in the Oort cloud with ai„ = 5 X 10 4 AU, ei„ 
0, ii n = 70° under the influence of the solar gravity and the galactic tide for the simple (dotted 
line) and the standard (solid line) models. Local density p = 0.1 Mq pc -3 is considered for 
both models. 

model (dotted line) the number of returns is only 4 during the integration time 4.5 x 10 9 
years. The difference between the two current models are also shown in the evolution 
of the inclination i, longitude of the ascending node !? and longitude of pericenter uj 
in Fig. 2. 

Fig. 3 represents the case when the initial inclination of the comet is %i n = 70°. 
Again, semi-major axis and eccentricity are dj„ = 5 x 10 4 AU, ej„ « 0. In this special 
case < di/dt > = for the simple model (see Eq. 3). For the standard model the orbit 
of the comet can be prograde, but also retrograde. The orbit of the comet is always 
prograde for !;„ < 90° . The orbit of the comet is always retrograde for ij„ > 90° . 

The number of returns of the comet to the inner part of the Solar System is 6, if 
ii n — 90°, both for the simple and standard models. 

For the difference between the number of cometary returns for the current models 
and more physical models compare Figs. 2, 3 with Figs. 5-8. 



4 Improved models 

Current models discussed in Sec. 3 must be treated as a rough approximation to reality. 
At present, an improved and more realistic physical access to galactic tide is in disposal 
(Klacka 2009a, 2009b). This section presents orbital evolution for the improved access 
to galactic tide. The results are compared to the evolution for the standard model. 
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4.1 Model I 

We are interested in the motion of the comet with respect to the Sun. The Sun is 
moving in a distance Rq — 8 kpc from the center of the Galaxy. Currently, the Sun 
is situated 30 pc above the galactic equatorial plane (Zq — 30 pc). Besides rotational 
motion with the speed (A — B) Rq the Sun moves with the speed 7.3 km/s in the 
direction normal to the galactic plane. Positional vector of the comet with respect to 
the Sun is r = (£, r), Q. 

Equation of motion is taken in the form (Klacka 2009, Eq. 29) 

d 2 i _ gm & ^ + B - ) [ A + B + 2j4cos (2 WQt )] £ 



di 2 



- 2A(A- B)sin(2 w t) r? 



+ (A - Bf (a/(6 2 + Zl) 1/2 + T 2 ) Rq Z a cos (w t) C , 

^Z = -^ST r l - 2A(A-B)Bin(2wo«) £ 
+ (A — B)[A + B — 2A cos (2 w t)] r) 
-{A- Bf (a/(6 2 + Zl) 1 ' 2 + T 2 ) R Zo sin (wot) C , 

|l = -^C- [4. G , + 2(^-5 2 )]C 
— 4 7r G £)' Zo [cos (wot) ^ — sin (wot) rf\ , 

d ^ = - [4nG Q + 2(A 2 B 2 )} Zo, 

r = + V 2 + C 2 , 

w = A - B , (6) 

where G is the gravitational constant, Mq is the mass of the Sun and the numerical 
values of the other relevant quantities are 

A = 14.25 km s _1 kpc -1 , 

B = - 13.89 km s _1 kpc -1 , 
A = 0.084 kpc" 1 , 
T 2 = 0.008 kpc~ 2 , 
g = 0.143 M pc~ 3 , 

Q =- 0.0425 M pc~ 3 kpc -1 , 

b = 0.25 kpc . (7) 

Eqs. (6) - (7) correspond to the model of the Galaxy presented by Dauphole et al. 
(1996). 

Fig. 4 shows the position of the Sun with respect to the plane of the galactic equator 
as a function of time for two improved models. The period of oscillations of the Sun 
with respect to the plane of galactic equator is 6.96 x 10 7 years for the Model I and 
7.26 x 10 7 years for the Model II. Maximal distance of the Sun is 87.8 pc from the 
plane of the galactic equator for the Model I and 91.14 pc for the Model II. 

The Sun is always located in the plane of galactic equator for the current models 
discussed in Sec. 3. 
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Fig. 4 Position of the Sun with respect to the plane of the galactic equator as a function of 
time. Galactic tides for the Model I (dotted line) and the Model II (solid line) are considered. 

Figs. 5 and 6 depict the orbital evolution of the comet as a numerical solution of 
Eq. (6), if also Eqs. (1) and (7) are taken into account. Comet is initially situated 
on the almost circular orbit (e^ n ~ 0) in the Oort cloud, = 5 x 10 4 AU. Orbital 
evolutions for two different initial inclinations are presented, ii n = 90° in Fig. 5 and 
ii n — 70° in Fig. 6. For comparison, the initial conditions of the comet are the same 
as for numerical calculations presented in Sec. 3. 

Figs. 5 and 6 represent the evolution of eccentricity e, inclination i, , longitude of 
the ascending node O, longitude of pericenter uj, perihelion distance q and aphelion 
distance Q during 4.5 x 10 9 years. Semimajor axes are practically constant during the 
whole integration time for both presented inclinations. 

Evolutions of the orbital elements in Figs. 5, 6 show the main differences between 
the current models described in Sec 3 and our improved Model I (compare with Figs. 
2, 3 ). 





Standard model 


Model I 




< q > [AU] 


< Q > [AU] 


< q > [AU] 


< Q > [AU] 


70° 


33792 


66889 


29484 


71526 


90° 


29827 


70685 


27574 


73605 


120° 


36614 


63922 


32934 


67884 



Table 1 Mean values of the perihelion and aphelion distances during the integration time 
4.5 X 10 9 years for the standard model and the Model I. Results for three different initial 
inclinations of the comet with a{ n = 5 x 10 4 AU, e ln « arc presented. 



Table 1 presents mean values of the perihelion and aphelion distances for the stan- 
dard model and the Model I. Various initial inclinations are presented. Both the stan- 
dard model and the Model I have minima of the mean values of the perihelion distances 
for the initial inclination ij n = 90° . Also the maxima of the mean values of the aphelion 
distances belong to the initial inclination ij„ = 90°. On the basis of our calculations 
we can state that the Model I shows lower mean values of the perihelion distances than 
the standard model and this holds irrespective to initial inclinations. The Model I also 
shows higher mean values of the aphelion distances than the standard model, again 
irrespective of the initial inclinations. 
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Fig. 5 Orbital evolution of the comet situated in the Oort cloud with cii n = 5 X 10 4 AU, 
ei n 0, ii n = 90° under the influence of the solar gravity and the galactic tide for the Model 
I. 



4.2 Model II 



Equation of motion is taken in the form (Klacka 2009a, Eqs. 26-27) 
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where G is the gravitational constant, M© is the mass of the Sun and the numerical 
values of the other relevant quantities are 

A = 14.2 km s _1 kpc -1 , 

B = — 12.4 km s _1 kpc^ 1 , 
A = 0.124 kpc~ 2 , 
T 2 = 1.586 kpc~ 4 , 
q = 0.130 M© pc~ 3 , 

q = - 0.037 M pc~ 3 kpc -1 , (9) 

see Eqs. (20)-(21) in Klacka (2009). If one wants to use other values of the Oort 
constants A and B, then he can use the following equation for calculation of mass 
density in the neighborhood of the Sun: 



Q = Qdisk + Qhalo , 

Bdlsk = 0.126 M© pc~ 3 , 

Qhalo = (4^G)- 1 [X {Galaxy) - X(disk) - X(bulge)] , 
X(Galaxy) = -(A- B) x (A + 35) 
X(disk) = 396.898 km 2 s~ 2 kpc~ 2 , 

X(bulge) = 0.625 km 2 s~ 2 kpc~ 2 . (10) 

The value X(Galaxy) = 611.800 km 2 s~ 2 kpc~ 2 holds for A = 14.2 km s _1 kpc -1 
and B = - 12.4 km s _1 kpc -1 . Moreover, Eq. (17) of Klacka (2009) can be used. 



Figs. 7 and 8 represent the orbital evolution of a comet as a numerical solution of 
Eq. (8), if also Eqs. (1) and (9) are taken into account. The comet is initially situated 
in the Oort cloud ( ai n = 5 x 10 4 AU) on the initially almost circular orbit as well as 
in pervious cases. Orbital evolutions for two different initial inclinations are presented, 
i in = 90° in Fig. 7, i in = 70° in Fig. 8. 

Figs. 7 and 8 depict the evolution of eccentricity e, inclination i, , longitude of the 
ascending node Q, longitude of pericenter u>, perihelion distance q and aphelion distance 
Q during the current lifetime of the Solar System. Semimajor axes are practically 
constant during the whole integration time for both presented inclinations. 

The evolutions of the orbital elements in Figs. 7 and 8 show the main differences 
between the current models described in Sec. 3 and our improved physical Model II 
(compare with Figs. 2 and 3 ). 





Standard model 


Model II 




< q > [AU] 


< Q > [AU] 


< q > [AU] 


< Q > [AU] 


70° 


33792 


66889 


25157 


75766 


90° 


29827 


70685 


27306 


73757 


120° 


36614 


63922 


33768 


66971 



Table 2 Mean values of the perihelion and aphelion distances during the integration time 
4.5 X 10 9 years for the standard model and the Model II. Three different initial inclinations of 
the comet with a,i n = 5 X 10 4 AU, e^„ « are considered. 
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Fig. 7 Orbital evolution of the comet with initial values aj n = 5 X 10 4 AU, ei n 0, ij. 
90° under the influence of the solar gravity and the galactic tide for the Model II 



Table 2 contains mean values of the perihelion and aphelion distances for the stan- 
dard model and the Model II. Various initial inclinations are presented as well as in 
Table 1. For the standard model, the minimum of the mean value of the perihelion dis- 
tance belongs to the initial inclination ij„ = 90°, but for the Model II the least value 
belongs to i in — 70° . Maximum of the mean value of the aphelion distance belongs 
to the initial inclination ij„ = 90° for the standard model and to ij„ = 70° for the 
Model II. The Model II shows lower mean values of the perihelion distances than the 
standard model for all of the presented initial inclinations. The Model II also shows 
higher mean values of the aphelion distances than the standard model for all of the 
presented initial inclinations. 
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Fig. 8 Orbital evolution of the comet with a^ n = 5 x 10 4 AU, e in si 0, i m = 70° under the 
influence of the solar gravity and the galactic tide for the Model II. 



5 Discussion 



Models I and II are based on more realistic physics than the standard or simple models. 
The Model I corresponds to the model of the Galaxy by Dauphole et al. (1986). This 
model of the Galaxy does not respect several important observational facts. It does 
not yield a flat rotation curve for great galactocentric distances. Mass density in the 
vicinity of the Sun is relatively large in comparison to newer values. Moreover, it does 
not produce the observed values of the Oort constants. Both of these observational 
facts are taken into account in the model of the Galaxy which served for our Model II 
(Klacka 2009a). This is the reason why we consider the Model II to be more relevant 
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than the Model I. In any case, both of these models take into account motion of the 
Sun in a more realistic way than it is considered in the simple and standard models. 

As a consequence of the previous discussion we can stress several facts. The stan- 
dard model yields mean value of the perihelion distance in w 35% higher than the 
real value for i in = 70°. This follows from Tables 1 and 2: < q > (standard model) 
/ < q > (Model I) = 1.146, < q > (standard model) / < q > (Model II) = 1.343. 
Moreover, the Model II is the only of the discussed models which yields (q)(ii n = 70°) 
< (q)(i ln = 90°). 

Let T q is a timescale on which a comet's perihelion changes in Aq. Let us consider 
Fig. 7 and Table 3 for finding r q for i k 2.3 x 10 9 years. The case of maximum in 
q — q(t) yields: 

r q (q = 4.8 x 10 4 AU, Aq = 0.4 x 10 4 AU; reality) = 9.4 x 10 7 yrs, 

while the relation presented by Levison and Dones (2007, p. 583) yields 

T q (q = 4.8 x 10 4 AU, Aq = 0.4 x 10 4 AU; literature) = 4.8 x 10 6 yrs. 

The mathematical relation yields 20-times smaller value than the real value. Now, let 

us consider minimum in q = q(t) for t ~ 2.3 x 10 9 years. The real value is 

T q (q = 1/4 x 10 4 AU, Aq = 1/2 x 10 4 AU; reality) = 4.7 x 10 7 yrs, 

while the mathematical relation yields 

T q (q = 1/4 x 10 4 AU, Aq = 1/2 x 10 4 AU; literature) = 2.6 x 10 7 yrs, 
or, using the value q — 38 AU from Table 3, 

Tq (q = 38 AU, Aq = 1/2 x 10 4 AU; literature) = 2.1 x 10 s yrs. 

Of course, smaller value of Aq may be taken into account. We can say that 

0.5 < Tq (literature) / T q (reality) < 20. 




1 1 1 years I 

Fig. 9 Comparison of the evolution of the mean perihelion distance of the comet with a; n = 
5 X 10 4 AU, e; n Ri 0, ij n = 90° for the standard model (dotted line), the Model I (dashed line) 
and the Model II (solid line). 



Table 3 presents minimal values of perihelion distances for the calculated models. 
While the standard model yields three cases with Qrnin 

< 50 AU, the Model I shows 
that the same situation occurs in one case, only. As we have discussed above, the 
relevant Model II yields six cases with q m i n < 50 AU. Table 4 presents values of the 
same quantities for initial galactocentric inclination 70 degrees. As a result we can 
state that the smallest values of q m i n exist for initial inclination corresponding to 90 
degrees. Fig. 9 depicts time average of evolution of perihelion distance (q) = (q(t)) = 
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Standard model 


Model I 


Model II 


q [AU] 


Q [AU] 


q [AU] 


Q [AU] 


q [AU] 


Q [AU] 


113.0 


101142.6 


1.5 


102594.9 


24.7 


102713.9 


38.7 


102137.5 


143.8 


102692.2 


101.3 


101620.5 


3.7 


101971.8 


115.7 


101718.5 


93.4 


102701.1 


71.9 


101056.8 


50.5 


102034.1 


104.2 


101700.1 


21.2 


101258.6 


187.5 


102263.4 


38.1 


102769.8 


111.6 


100130.6 


194.9 


101536.0 


25.9 


102648.5 






192.1 


101758.7 


43.7 


102744.2 






90.5 


101707.7 


1.0 


102455.4 






149.0 


102142.4 


37.5 


102046.6 








101552.9 







Table 3 Extremal perihelion and aphelion distances of the comet with a^ n — 5 x 10 4 AU, 
ei n « 0, ii n = 90° during the integration time 4.5 X 10 9 years for the standard model, the 
Model I and the Model II. 



Standard model 


Model I 


Model II 


q [AU] 


Q [AU] 


q [AU] 


Q [AU] 


q [AU] 


Q [AU] 


3465.38 


98209.0 


2883.6 


99568.6 


2968.3 


99308.4 


3036.7 


98793.8 


4539.5 


97433.9 


5067.2 


96652.2 


3063.5 


98775.0 


4168.9 


98050.3 


3446.5 


98680.8 


4466.9 


96966.1 


3190.1 


99294.7 


5299.0 


96224.7 


3062.8 


98648.0 


3689.2 


98594.1 


3267.2 


98839.0 






2714.9 


100020.7 


5441.0 


95944.8 






5114.2 


96658.8 


4385.5 


97489.7 






2753.9 


99851.9 


3970.5 


97988.4 






3407.8 


99033.9 


5573.5 


95894.2 



Table 4 Extremal perihelion and aphelion distances of the comet with aj n = 5 X 10 4 AU, 
ei n ~ 0, ii n = 70° during the integration time 4.5 X 10 9 years for the standard model, the 
Model I and the Model II. 



J q(t')dt' I t. The evolution of perihelion distance yields minimal values for the best 
physical model represented by our Model II. 

The effect of Solar System's vertical motion above and below the galactic plane on 
the possible terrestrial mass extinctions was discussed in the past (see, e.g., Weissman 
1990, Rampino and Stothers 1984, Schwartz and James 1984). Results corresponding 
to Fig. 4 show that period of the Solar System's vertical oscillations is 72.6 x 10 6 
years. The last Sun's position in the galactic equatorial plane occured before 3.87 X 
10 6 years. The last extinctions had to occur before (3.87 x 10 6 —ex a 3 / 2 ) years, where 
semi- major axis a of the comet is given in astronomical units and < e < 1. The value 
of e describes the initial position of the comet on its orbit (e = 1/2 corresponds to 
aphelion). Since the Solar System's oscilations are periodic, other possibilities are [(n 
x 72.6 / 2 + 3.87 ) x 10 6 —ex a 3 / 2 ] years, where n is an integer number, which 
multiplies the half period of the Solar System's oscilations. Thus, also (40.2 x 10 6 — e 
x a 3/2 ) years (n = 1), or, (76.5 x 10 6 - e x a 3/2 ) years (n = 2) in the past. The values 
could correspond to the value of about 65 x 10 6 years of the late Cretaceous extinctions 
which included the disappearance of the dinosaurs (Rampino and Stothers 1984). The 
mass extictions on the Earth have periodic behaviour with the period approximately 
33 x 10 6 years. The case n = 2 could occur for a = 5.1 x 10 4 AU / e 2/3 . The comet 
was increasing its heliocentric distance for some part of its motion before the hit of 
the Earth, if 5.1 x 10 4 AU < a < 8.1 x 10 4 AU. If the comet was still moving toward 
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the Sun, then e < 1/2 and a > 8.1 x 10 4 AU. The case n = 1, for 33 x 10 6 years 
in the past, yields a = 3.73 x 10 4 AU / e 3/2 . If e < 1/2, then a > 5.9 x 10 4 AU. 
Similarly, the case n = 3, for 99 x 10 6 years in the past, yields a = 5.74 x 10 4 AU 
/ e 3 / 2 . Although the galactic tidal forces exclude stability of the Oort cloud for large 
values of semi-major axis, a > 8 x 10 4 AU, a close encounter of a star or an interstellar 
cloud could disturb a comet of the Oort cloud. Finally, there is also a possibility that 
objects of interstellar origin hit the Earth. 

There is an another feature of the effect of Solar System's vertical motion above 
and below the galactic plane on the possible terrestrial mass extinctions. We should 
await, on the basis of Fig. 7 and Table 3, that the period of possible exceptionally 
strong mass extinction is that the period of possible mass extinctions is (4.5-5.0) x 10 8 
years, if the galactic tide is the reason of the mass extinctions. This is twice the value 
presented by Schwartz and James (1984). It seems that we are not able to explain the 
major extinction events on the basis of the galactic tide alone. 

6 Conclusion 

The paper treats the effect of the galactic tide on a cometary motion with respect to 
the Sun. Our detailed numerical calculations show that the F— terms play an important 
role in cometary orbital evolution. The J 1 — terms cause that periods of secular evolution 
of the orbital elements differ from the periods of conventional models. The real period 
is 1.7-times smaller than the period of the conventional models, if the initial cometary 
inclination is close to 90 degrees, eccentricity is close to zero and semi-major axis is 
about 5 x 10 4 AU. Shorter periods of secular evolution of orbital elements lead to the 
conclusion that frequency of cometary returns into the inner part of the Solar System 
is about two times higher than the frequency following from the conventional models. 

A comet from the Oort cloud can be a source of the dinosaurs extinction at about 
65 Myr ago. Either a close encounter of a star or an interstellar cloud disturbed a comet 
in the Oort cloud in the way that its semi- major axis increased/decreased above the 
value 5.1 x 10 4 AU and the comet hit the Earth. Or, an object of interstellar origin 
hit the Earth. 
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